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ABSTRACT:  Amine,  NH2,  is  of  interest  as  a  linker  between 
organic  molecules  and  graphene  in  novel  biotechnologies.  We 
used  ab  initio  electronic  structure  calculations  to  study  NH2 
adsorption  on  graphene  in  the  presence  of  surface  defects  and 
other  adsorbates,  including  N,  B,  F,  H,  and  OH.  Amine  is 
found  to  form  a  semi-ionic  bond  of  0.778  eV  on  pristine 
graphene.  Its  binding  is  found  to  be  modified  near  other 
defects,  and  the  adsorption  energy  is  dependent  on  the 
neighbor  order. 


■  INTRODUCTION 

In  recent  years  graphene  has  become  one  of  the  most 
interesting  materials  in  nanotechnology  in  part  because  it 
opens  opportunities  in  a  range  of  technologies,  including 
sensors,  electronics,  and  hydrogen  storage.  However,  some  of 
the  properties  of  graphene  must  be  modified  for  it  to  be  useful 
in  these  applications.  Surface  chemical  functionalization  is 
particularly  effective  at  changing  the  properties  of  graphene 
because  it  is  entirely  made  of  surface  atoms.  A  wide  range  of 
surface  chemistries  are  being  investigated  to  chemically 
functionalize  graphene  and  to  alter  its  properties.1-4 

One  class  of  functionalizations  involve  amines,  the  simplest 
being  NH2.  Amines  are  of  interest  because  they  are  used  in 
catalytic  asymmetric  synthesis,5  chromatography,6  and  as  a 
biconjugate  linker7  (i.e.,  between  graphene  and  DNA).  Previous 
calculations  have  shown  that  NH2  binds  at  the  on-top  sites  of 
pristine  graphene  with  binding  energies  and  mobilities  that 
would  allow  it  to  diffuse  readily  across  the  surface.8  It  is  known 
that  adsorbates  and  surface  defects  can  modify  molecular 
adsorption  properties.9'10  Previous  experimental11  and  theoreti¬ 
cal  work12  on  carbon  nanotubes  and  graphene  have  shown  that 
molecular  adsorption  energies  often  are  increased  at  defects 
such  as  those  related  to  oxygen  and  that  these  defects  can 
modify  sensor  behaviors. 

Here  we  consider  NH2  adsorption  on  graphene  in  the 
presence  of  boron  or  nitrogen  substitutional  defects  as  well  as 
fluorine,  hydrogen,  hydroxyl,  or  other  amine  adsorbates.  We 
chose  to  look  at  these  defects  because  depending  on  the  growth 
processes  involved16  they  can  appear  in  amine— graphene 
systems.  We  discuss  the  computational  methods  used,  the 
effects  of  a  single  defect,  and  then  how  neighboring  defects 
affect  the  adsorption  properties  of  NH2  on  graphene. 

■  METHODS 

To  find  the  ground-state  structures  of  these  graphene  systems, 
we  used  Quantum  Espresso17'18  to  do  spin-unrestricted  plane- 
wave  density  functional  theory  (DFT)  calculations.  We  used 
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the  Perdew— Burke— Ernzerhof  (PBE)  exchange-correlation 
parametrization  of  the  generalized  gradient  approximation 
(GGA)  with  Vanderbilt  ultrasoft  pseudopotentials.  To  include 
weak  dispersive  interactions  (van  der  Waals  forces),  semi- 
empirical  dispersion  terms  (DFT-D)  were  included.19-21  12- 
by-12  rectangular  graphene  supercells  and  an  interlayer  spacing 
of  15  A  were  used. 

Figure  la  gives  the  relative  positions  of  two  graphene  defects 
in  our  models.  The  primary  defect  is  placed  near  the  middle  of 
the  graphene  sheet  at  position  0.  The  primary  defect  may  be  a 
substitutional  B  or  N  or  an  adsorbate:  F,  H,  OH,  or  NH2.  A 
second  defect,  always  NH2,  is  adsorbed  near  the  primary  defect 
at  one  of  the  sites  labeled  1—5.  The  positions  labeled  1,  2,  and  3 
are,  respectively,  known  as  the  ortho,  para,  and  meta  positions. 
In  this  work,  when  both  defects  are  adsorbates  they  are  both  on 
the  same  side  of  the  graphene  sheet. 

■  RESULTS  AND  DISCUSSION 

Primary  Defects.  We  first  consider  NH2  as  a  primary  defect 
adsorbed  alone  on  pristine  graphene.  We  find  that  it  adsorbs  at 
an  on-top  site,  as  shown  in  Figure  lb.  This  results  in  the  N 
atom  being  directly  over  the  C  atom,  to  which  it  is  bound  with 
a  C— N— H  angle  of  107.8°  and  with  each  H  atom  pointing 
toward  the  center  of  neighboring  benzyl  rings.  This  arrange¬ 
ment  of  atoms  is  consistent  with  the  polar  nature  of  the  amine 
molecule;  the  lone  pair  electrons  of  the  N  atom  are  oriented 
away  from  one  of  the  neighboring  pz  orbitals.  In  this 
configuration,  the  C— N  bond  length  is  1.518  A,  whereas  the 
C  atom  to  which  the  NH2  molecule  is  attached  is  0.583  A  out 
of  the  graphene  plane.  The  C— N  bond  length  indicates  that  the 
bond  is  semi-ionic.  We  find  that  the  NH2  has  an  adsorption 
energy  of  0.778  eV,  consistent  with  results  given  by  Milowska 
et  al.  The  reason  for  the  small  adsorption  energy  is  that  NH2 
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Figure  1.  (a)  Section  of  the  rectangular  graphene  supercell  where  the 
numbers  on  the  atoms  represent  the  neighbor  order  with  respect  to 
the  primary  defect  position  0  in  the  graphene  lattice.  In  this  example, 
the  primary  defect  is  a  substitutional  N  atom.  The  atoms  labeled  1,  2, 
and  3  are,  respectively,  in  the  ortho,  meta,  and  para  positions  in 
relation  to  position  0.  (b)  An  amine  group  is  shown  on  the  graphene 
surface  where  the  blue  atom  is  a  N  atom  and  the  white  atoms  are  H 
atoms. 


has  an  electron  that  it  would  like  to  donate,  whereas  the  carbon 
atoms  in  graphene  already  have  full  shells. 

To  understand  the  nature  of  the  bond  between  graphene  and 
NH2,  we  examined  their  electronic  density  of  states  (DOS). 
Figure  2a  gives  the  DOS  of  defect-free  graphene  and  the 
projected  density  of  states  (PDOS)  of  isolated  NH2. 
Comparing  the  relative  energies  of  the  highest  occupied 
molecular  orbital  (HOMO)  and  lowest  unoccupied  molecular 
orbital  (LUMO)  states  of  NH2  with  the  Fermi  energy  of 


DOS 


Figure  2.  Density  of  states  of  graphene— NH2  systems,  (a)  Spin- 
polarized  DOS  of  isolated  graphene  (red)  and  isolated  NH2  (black  and 
blue  (electron  spin  up  and  spin  down,  respectively))  at  infinite 
separation.  The  black  and  blue  labels  situated  to  the  right  of  the  peaks 
indicate  which  atomic  orbital (s)  contribute  most  to  a  particular  NH2 
molecular  orbital,  (b)  DOS  of  the  graphene  sheet  with  one  NH2 
molecule  adsorbed  on  it  is  shown  in  red.  PDOS  of  the  pz  orbital  of  the 
NH2  N  atom  is  blue.  PDOS  of  the  pz  orbital  of  the  C  atom  to  which 
the  N  atom  is  bounded  is  green. 


graphene,  we  can  determine  if  in  the  combined  system  charge 
should  be  transferred  via  band  alignment.23  If  the  HOMO  level 
of  NH2  is  above  the  Fermi  energy  of  graphene,  then  charge 
transfers  from  NH2  into  the  graphene  conduction  band.  If  the 
NH2  LUMO  is  below  the  graphene  Fermi  energy,  then  charge 
transfers  from  the  graphene  valence  band  into  the  NH2  LUMO 
states.  Our  results  show  that  the  HOMO  of  isolated  NH2  is  ~2 
eV  below  the  Fermi  energy  of  graphene,  whereas  the  NH2 
LUMO  is  several  electro nvolts  above  the  Fermi  energy.  Thus, 
in  the  NH2— graphene  system  any  charge  transfer  is  not  due  to 
band  alignment.  We  find  that  the  NH2  HOMO  is  located  on 
the  Np  orbital,  in  an  energy  region  where  the  carbon  pz  orbitals 
(in  the  form  of  /r-bonds)  dominate.  This  suggests  that  a 
covalent  C— N  bond  may  be  formed  between  graphene  and 
NH2.  Figure  2b  shows  the  effect  that  the  bonding  has  on  the 
PDOS.  When  the  NH2  molecule  is  adsorbed  on  graphene  a 
defect  peak  appears  at  the  Fermi  level,  disrupting  the  Dirac 
point.  The  PDOS  shows  that  it  is  due  to  a  change  in  the  pz 
orbital  of  the  C  atom  to  which  the  adsorbate  is  bound  with  a 
corresponding  peak  from  the  Npz  orbital.  Thus,  whereas  NH2 
would  rather  donate  an  electron  and  the  C  atoms  have  full 
shells  via  the  o  and  n  bonds,  a  C  atom  and  an  NH2  molecule 
can  form  a  weak  bond  between  the  Np  and  Cp  .  These  features 
are  consistent  with  the  C— N  bond  being  semi-ionic. 

We  now  consider  the  other  primary  defects  of  interest  here. 
Table  1  gives  the  binding  energies  of  each  primary  defect,  the 


Table  1.  Adsorption  Energy,  Charge  Transfer,  and 
Adsorbate  Bond  Length  for  the  Primary  Defects** 


defect 

energy 

[eV] 

adsorbate  charge 

[e] 

C  charge  [e] 

bond  length  [A] 

B 

5.799 

0.115 

NA 

NA 

N 

7.017 

0.011 

NA 

NA 

F 

2.310 

-0.394 

0.228 

1.678 

H 

1.324 

0.224 

-0.158 

1.110 

nh2 

0.778 

0.030 

0.108 

1.518 

OH 

1.227 

-0.136 

-0.187 

1.407 

"Positive  adsorbate  charge  value  means  electrons  have  been  donated 
by  the  atom  (molecule),  leaving  a  positively  charged  entity.  B  and  N 
are  substitutional  defects  in  the  graphene  lattice  and  thus  do  not  have 
C  atom  charge  difference  or  bond  lengths  from  C.  The  bond  lengths 
reported  for  NH2  and  OH  are  the  C— N  and  C— O  distances, 
respectively. 


Lowdin  charge  transfer  into  or  out  of  the  defect,  and  the  C— X, 
where  X  =  {F,  H,  O  (in  OH),  N  (in  NH2)},  the  bond  length 
between  an  adsorbate  and  the  C  atom  it  is  attached  to.  For  a 
substitutional  defect  the  formation  energies  in  Table  1  were 
computed  by 

A  E  =  (£gry  +  Ec )  —  ( £GR  +  Ey)  (i) 

where  EGR  is  the  energy  of  the  graphene  supercell,  EY  is  the 
energy  of  a  lone  N  (or  B)  atom,  £GR  Y  is  the  energy  of  the 
graphene  supercell  with  the  N  (or  B)  atom  replacing  one  of  the 
C  atoms,  and  Ec  is  the  energy  of  a  lone  C  atom.  A  negative  AE 
indicates  an  exothermic  process.  In  cases  when  the  primary 
defect  is  an  adsorbate,  the  energy  was  computed  using 

AE  =  Egrz  —  (Egr  +  Ez)  (2) 

where  EGR  z  is  the  energy  of  the  combined  graphene  and 
adsorbate  system  and  Ez  is  the  energy  of  the  adsorbate.  Among 
the  adsorbates  studied  here,  F  has  the  largest  adsorption  energy 
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when  adsorbed  alone  and  NH2  has  the  smallest.  The  B  and  N 
substitutional  defects  are  found  to  be  exothermic,  consistent 
with  expectations.24  Comparing  the  C— X  bond  lengths  of  the 
primary  adsorbate  defects  with  the  accepted  covalent  bond 
lengths,25  we  find  that  all  of  them  are  consistent  with  the 
bonding  being  semi-ionic. 

In  Figure  3,  we  illustrate  the  effect  that  adsorbate  defects 
have  on  the  calculated  charge  of  the  C  atoms  in  the  graphene 


Figure  3.  Effect  of  adsorbates  on  the  charge  in  surrounding  C  atoms. 
The  neutral  value  is  set  to  the  (Lowdin)  charge  on  the  C  atoms  in 
clean  graphene,  and  differences  from  that  value  are  colored  in  shades 
of  red  or  blue.  The  dark-gray  atoms  show  the  position  of  adsorbates. 
The  color  of  the  C  atom  to  which  the  adsorbate  is  attached  is  not 
indicative  of  the  charge  difference  on  that  atom;  when  F  is  attached, 
the  C  atom  has  a  charge  difference  of  0.23  e,  when  NH2  is  adsorbed, 
the  charge  difference  on  the  C  atom  is  0.11  e;  all  other  C  atoms  are 
represented  by  the  color  on  the  atom,  (a)  F  adsorption,  (b)  OH 
adsorption,  and  (c)  NH2  adsorption. 


sheet.  There  the  colors  of  the  C  atoms  represent  the  change  in 
the  charge.  Charge  transfer  is  due  to  the  difference  in 
electronegativities,  which  gives  a  charge  transfer  between  the 
C  atom  and  the  adsorbate.  In  Figure  3a  for  F  adsorption,  we  see 
that  more  than  just  the  C  atom  bonded  to  the  F  atom  is 
involved;  while  most  of  charge  difference  is  in  that  C  atom,  the 
nearby  C  atoms  also  show  small  charge  transfers.  The  charge 
transfer  on  the  neighboring  C  atoms  is  an  order  of  magnitude 


smaller  than  that  on  the  C  to  which  the  F  is  bound.  From 
columns  3  and  4  of  Table  1,  we  see  that  the  charge  removed 
from  C  atom  is  less  than  the  charge  transferred  to  the  F  atom, 
and  thus  charge  must  be  pulled  from  other  nearby  carbon 
atoms. 

The  charge  transfer  and  polarization  in  Figure  3  may 
influence  clustering  of  multiple  adsorbates  and  contribute  to 
other  effects,  such  as  distortion  of  the  lattice.26  Figure  3a  shows 
a  triangular  symmetry  of  the  charge  distribution  within  the 
graphene  sheet  around  the  F  adatom.  Triangular  symmetry  has 
been  related  to  the  low  energy  of  clustering  of  H  on  a  graphene 
sheet  by  Casolo  et  al.27  It  is  likely  that  the  lack  of  triangular 
symmetry  seen  in  Figure  3b, c  for  OH  and  NH2  adsorption 
suggests  that  these  molecules  may  exhibit  different  high 
concentration  adsorption  conformations  than  those  for  H  and 
F. 

Secondary  Defects.  Figure  4  gives  the  adsorption  energies 
of  a  secondary  NH2  in  the  presence  of  other  (primary) 


Sublattice 

A  B  A  B  B  A 


Neighbor  Order 

Figure  4.  Adsorption  energies  of  the  secondary  NH2  as  a  function  of 
neighbor  order  from  the  primary  defects  listed  in  the  Figure.  The 
dashed  line  represents  the  adsorption  energy  NH2  molecule.  The  gray 
region  shows  an  analytical  description  based  on  the  nearest-neighbor 
tight  binding  model  of  graphene,  as  described  in  the  text.  The  axis 
labeled  “Sublattice”  indicates  the  sublattice  (A  or  B)  of  the  C  atom  to 
which  the  secondary  defect  is  attached. 


defects.28  The  energies  are  given  by  eq  2,  where  £GR  now 
represents  the  energy  of  the  combined  graphene  and  primary 
defect  system  and  EGR  z  is  the  energy  of  the  combined 
graphene,  primary  defect,  and  NH2  system.  Thus,  the  total 
energy  of  the  primary  defect  adsorbed  on  graphene  is  used  as 
the  zero  of  the  energy  when  a  secondary  NH2  molecule  is 
adsorbed.  The  secondary  NH2  adsorption  energies  have  similar 
dependencies  on  separation  from  the  several  primary  defects. 
The  separation  dependence  is  shaped  primarily  by  the 
electronic  structure  of  graphene's  conjugated  n  electrons,  as 
we  discuss  below.  The  relative  energies  when  both  adsorbates 
are  on  the  same  side  of  the  graphene  sheet  are  similar  to  the 
trend  that  Casolo  et  al.  give  in  table  2  of  ref  27  for  two  H  atoms 
adsorbed  on  the  same  side  of  the  sheet. 

Substitutional  nitrogen  defects  may  occur  whenever  amines 
are  adsorbed  onto  graphene.  NH2  does  not  form  a  stable  bond 
directly  with  the  substitutional  N.  This  is  due  to  both  N  atoms 
in  the  bond  wanting  to  donate  an  electron.  The  adsorption 
energy  of  the  NH2  on  a  first-  or  third-order  neighbor  site  of  the 
substitutional  N  is  nearly  37%  larger  than  the  adsorption  energy 
of  a  singly  adsorbed  NH2.  The  adsorption  energy  of  the 
secondary  NH2  when  it  is  on  the  second-  or  third-order 
neighbor  sites  is  closer  to  the  singly  absorbed  NH2  energy.  The 
secondary  adsorption  energy  of  an  NH2  in  the  presence  of 
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another  NH2  is  consistently  lower  than  when  NH2  adsorbs  near 
a  substitutional  N.  This  suggests  that  in  experiment  we  should 
expect  to  see  more  NH2  adsorbed  onto  the  graphene  when 
substitutional  N  are  present  than  when  only  NH2  adsorbates 
are  present. 

To  understand  the  variation  of  adsorption  energy  of 
secondary  NH2  with  distance,  we  recall  that  the  primitive  cell 
of  graphene  contains  two  carbon  atoms  and  forms  a  bipartite 
lattice.29’30  Lattice  sites  related  to  these  two  different  atoms  of 
the  primitive  cell  by  symmetry  operations  form  two  sublattices, 
A  and  B.  From  DFT  results  in  Figure  4,  we  notice  that 
adsorption  energy  changes  depending  on  whether  the  defects 
are  on  the  same  sublattice,  for  example,  both  on  A  sublattice 
sites,  or  are  on  different  sublattices.  The  DFT  calculations 
include  electron  orbitals  in  a  wide  range  of  energies  including 
many  effects,  for  example,  elastic  deformation  of  the  lattice. 
However,  only  orbitals  in  a  narrow  range  of  energies  near  the 
chemical  potential  that  form  the  conduction  network 
(conjugated  n  electrons)  are  substantially  affected30  by  the 
bipartite  nature  of  the  lattice.  This  suggests  that  conduction 
electrons  play  the  dominant  role  in  defining  the  variation  of 
adsorption  energy  with  distance  between  the  two  defects. 

To  demonstrate  this,  we  model  the  conjugated  n  electrons  by 
a  tight-binding  Hamiltonian,  in  which  electrons  from  carbon  pz 
orbitals  hop  between  neighboring  pz  orbitals  of  carbon  with 
amplitude  t.29  Adsorbates  disrupt  this  hopping  by  altering  local 
electronic  configurations,  thus  introducing  scattering  centers. 
The  simplest  distance-dependent  contribution  to  the  adsorp¬ 
tion  energy  of  a  pair  of  adsorbates  involves  two  such  scattering 
events  (vertices)  and  is  given  by  a  polarization  bubble  n(R), 
where  R  is  the  distance  between  the  two  vertices.  Because  of 
bipartite  nature30  of  the  graphene's  lattice,  n(R)  depends  on 
the  sublattice  index 


dco  dq  dq  /R(q-q/) 

(2/r)5 


G^ico,  q )GJi(ia),  q') 


(3) 


where  i  and  j  refer  to  A  or  B  sublattices.29  Near  the  half- filling, 
the  spectrum  of  graphene  has  two  nonequivalent  conical  valleys 
K  and  K'.30  The  characteristic  momentum  difference  K  —  K' 
between  the  two  adjacent  valleys  gives  rise  to  fast  spatial 
oscillations.  The  momentum  integration  in  eq  3  runs  over  the 
entire  Brillouin  zone;  however,  at  (or  near)  the  half-filling  the 
dominant  contributions  come  from  the  regions  near  the  conical 
points  K(K').  The  lattice  Green's  functions  on  different 
sublattices  Gab(  ico,  K  +  q)  and  GAB(iu>,  K'  +  q)  acquire  an 
additional  phase  factors  of  ( qx  +  iqy)/q  and  ( qx  —  iqy)/q , 
respectively,  when  Iql  IKI 29,31  This  sublattice-dependent 
phase  factors  together  with  e*R^q-q/')  produce  an  overall 
oscillatory  factor.  In  particular,  when  the  adsorbates  are  on 
the  same  sublattice 


eiR(q  q')  |e»RK  +  e<RK'|2  =  2  +  2  COS  R(K  -  K')  (4) 

whereas  when  the  adsorbates  are  on  the  different  sublattices 

q  ±  iq 

iR(q-q')  x  /  |  i6+in/2  iRK  .  -iO-in/2  iRK  ,2 

e  vn  n  / -  — »  le  e  +  e  el 

=  2  +  2  cos{R(K  -  K')  4-  20  +  n)  (5) 

where  6  is  the  angle  between  R  and  K  —  K.  In  addition,  the 
particle-hole  symmetry  of  the  Hamiltonian  at  half-filling 
requires  and  TI^  to  have  opposite  signs,  as  discussed  in 
detail  by  Saremi.31  The  magnitude  of  eq  3  depends  on  details, 


such  as  strength  of  the  scattering,  polarization  of  the  bonds,  and 
chemical  potential,  and  is  cumbersome  to  evaluate.  However, 
the  above  symmetry  relation  together  with  eqs  4  and  5  provides 
a  qualitative  description  of  the  oscillatory  behavior  found  with 
DFT  calculations.  In  Figure  4,  the  overall  magnitude  of  the 
analytical  result  is  adjusted  to  match  the  DFT  data.  Effects  due 
to  lower-laying  orbitals,  such  as,  for  example,  local  strain  in 
graphene  lattice,  is  likely  to  play  a  role  in  deviation  of  the 
adsorption  energies  from  the  energy  given  by  the  above 
qualitative  description.  This  can  be  especially  pronounced  at 
short  distances,  where  local  deformation  of  C— C  bonds  can 
play  an  important  role.26 

The  C-X,  where  X  =  {F,  H,  O  (in  OH),  N  (in  NH2)},  bond 
lengths  for  primary  defects  and  the  secondary  NH2  for  different 
distances  between  defect  sites  are  shown  in  Figure  5.  In  Figure 
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A  B  A  B  B  A 


Figure  5.  Bond  lengths:  (a)  The  bond  length  between  the  primary 
adsorbate  (C—  bond)  and  carbon  atom  to  which  it  is  bound.  This 
excludes  the  substitutional  defects,  (b)  The  bond  length  between 
carbon  and  the  N  atom  in  the  secondary  NH2  (C— N  bond).  NH2  can 
be  both  a  primary  and  secondary  defect.  The  Nsub  and  Bsub 
designations  are  for  the  substitutional  N  and  B  atoms  and  thus  do  not 
have  values  in  panel  a. 


5a,  we  give  the  values  for  the  bond  lengths  for  the  primary 
defects  adsorbed  on  top  of  a  carbon  atom.  The  bond  lengths 
for  the  secondary  NH2  (for  the  bond  between  nitrogen  and 
carbon)  are  given  in  Figure  5b.  The  variation  in  the  C— N  bond 
length  for  the  secondary  NH2  defect  placed  at  different 
positions  relative  to  the  primary  defect  (the  first  through  the 
fifth  order  neighbor  sites)  is  <5%.  The  C— X  bond  changes  as 
much  as  25%  with  distance  from  the  NH2  molecule.  This 
variation,  together  with  the  variation  in  charge  transfer 
discussed  below  affects  the  scattering  of  conduction  electrons 
on  primary  and  secondary  defects  and  thus  can  contribute  to 
the  variation  of  the  adsorption  energies  of  the  secondary  NH2. 

Figure  6a, b  show  the  total  (Lowdin)  charge  transferred  into 
the  primary  and  secondary  defects,  respectively,  when  two 
defects  are  present.  The  charge  transferred  into  the  primary 
defect,  of  any  particular  type  should  be  compared  with  the 
single  defect  value  listed  in  Table  1  for  the  same  type  of  defect. 
Similarly,  the  charge  transferred  into  the  secondary  defect, 
always  an  NH2  molecule,  should  be  compared  with  the  single 
NH2  value  listed  in  Table  1  for  NH2.  The  values  shown  in 
Figure  6  have  large  deviations  from  the  single  defect  cases.  The 
average  deviation  of  the  charge  transfer  is  about  44  and  80%  for 
the  primary  and  secondary  defects,  respectively.  The  largest 
deviations  are  observed  in  the  cases  where  the  secondary  NH2 
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Figure  6.  Lowdin  charge  transfer  into  each  defect  when  NH2  is 
adsorbed  near  a  primary  defect,  (a)  Charge  transfer  into  the  primary 
defect,  (b)  Charge  transfer  into  the  NH2  defect.  A  positive  value 
means  electrons  have  been  donated  by  the  atom  (molecule);  leaving  a 
positively  charged  entity. 


is  adsorbed  on  a  substitutional  defect.  In  these  cases  the 
deviation  can  be  attributed  to  the  fact  that  NH2  forms  a  direct 
bond  with  substitutional  defect  instead  of  a  carbon.  There  is 
also  large  deviation  in  the  charge  transferred  when  NH2  is 
located  near  (first  or  second  neighbor)  an  OH  or  NH2.  In  the 
case  of  two  adsorbed  NH2  groups;  or  a  NH2  and  OH  pairing; 
the  large  deviation  in  charge  transfer  is  likely  partially  due  to 
inaccuracies  in  computing  the  Lowdin  charge.  In  these  cases; 
the  adsorbed  molecules  are  rotated  around  their  C-X  bonds  in 
such  a  way  as  to  indicate  hydrogen  bonding.  Hydrogen  bonding 
can  affect  the  polarization  of  the  adsorbate  group;  thus 
changing  the  distribution  of  charge  within  the  group.  In  this 
case  the  Lowdin  calculation  does  not  attribute  some  of  the 
charge  to  the  atoms;  where  it  is  expected  to  reside.  Overall;  we 
find  that  the  difference  in  charge  transfer  for  the  cases  of  a  pair 
of  defects  and  a  lone  defect  decreases  as  neighbor  order 
increases  and  defects  become  more  isolated  from  each  other. 

■  SUMMARY 

In  summary;  we  studied  the  bonding  of  NH2  molecules  near 
other  graphene  defects.  These  defects  include  adatoms; 
admoleculeS;  and  substitutional  defects.  The  binding  energy 
of  an  NH2  depends  on  its  distance  from  a  primary  defect.  The 
distance  dependence  is  due  to  scattering  processes  in 
graphene's  conjugated  n  electrons.  At  a  particular  distance 
from  the  primary  defect  the  adsorption  energy  of  NH2  also 
depends  on  the  type  of  primary  defect.  We  find  that  NH2 
chemisorbs  to  the  graphene  surface  with  a  semi-ionic  bond. 
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